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Abstract 

We consider the problem of spectra for acoustic type of turbulence gen- 
erated by shocks being randomly distributed in space. We show that for 
turbulence with a weak anisotropy such spectra have the same dependence 
in A;-space as the Kadomtsev-Petviashvili (KP) spectrum: E{k) ~ k^"^. How- 
ever, the frequency spectrum has always the falling ~ independently 
on anisotropy. In the strong anisotropic case the energy distribution relative 
to wave vectors takes anisotropic dependence forming in the large k region 
the spectra of the jet type. 

PACS: 52.35.Ra, 52.35.Tc, 47.35.Rs, 47.35.Jk 

1 Introduction 



The acoustic type turbulence is ubiquitous in space and laboratory plasmas. Typ- 
ical example represents MHD turbulence (for review see [1] and references therein) 
in the presence of external magnetic field at a moderate level of (3 {(3 is the ra- 
tio of kinetic plasma pressure to the magnetic pressure). When wave amplitudes 
are small the turbulence can be described as an ensemble of linear noninteracting 
waves with their frequencies Uk and wave-vectors k. In the long-wave limit wave 
frequency can be considered as linear function of a wave number: ujk ~ k. In 
the thermodynamic equilibrium limit these waves are distributed according to the 
Jeans law so that the energy of each wave = T where T stands for the "wave's" 
temperature. When the wave amplitudes increase (with the growth of the energy 
source) we first come to the regime of weak turbulence. The nonlinear effects are 
weak in comparison with linear wave dispersion but statistical characteristics of 
the system change significantly. In the weak turbulence regime, in addition to the 
thermodynamic distribution solution, it emerges the additional solution having the 
spectrum of the Kolmogorov type: 

e{k) = C {pPcsf^^ k-'/^ 



1 



where P is the energy flux toward small scales, p the density, the sound ve- 
locity and C is the dimensionless constant of the order of one. In the spherical 
normalization this energy distribution reads as follows 



This spectrum was found in 1971 by Zakharov and Sagdeev [2] as the exact solution 
of the wave kinetic equation for acoustic waves in isotropic medium, when in the 
long- wave region the dispersion relation is linear Uk ^ Cgk (for details, see [3]). In 
this case the criterion of weak turbulence that is determined as the weakness of 
the nonlinearity with respect to the wave dispersion, is written as 



where Ap is the fluctuating part of the density p, and A is the dispersion length 
0. For larger amplitudes when the nonlinear effects become comparable or larger 
the wave dispersion this criterion ([2]) breaks down. As it is well known in gas dy- 
namics, in this case the main nonlinear effect for acoustic waves is nothing else but 
the wave breaking that results in the formation of shocks (or discontinuities). It is 
well known also that this process in compressible flows can be treated in terms of 
the formation of folds in the classical catastrophe theory [5]. In the gas-dynamic 
case, breaking areas can be completely characterized using the mapping defined by 
the transition from the usual Eulerian to the Lagrangian description. A vanishing 
of the Jacobian of the mapping corresponds to the emergence of a singularity for 
the spatial derivatives of the velocity and density of the fluid. Physical meaning 
of this effect corresponds to the intersection of particle trajectories. In the general 
situation first time the Jacobian vanishing happens in one isolated point. In colli- 
sionless plasma this process can continue and leads further to form the multi-flow 
region expanding in the transverse (relative to the main flow) direction according 
to the following scaling law R± ~ y/to — t, where to stands for the first moment 
when Jacobian J turns to zero, while the region width (in the longitudinal direc- 
tion, along the main flow) increases more slowly, R\\ ~ (to — t)^''^ (see, e.g. [HI [7]). 
Thus, the result of the breaking consists in the formation of structures in the form 
of pancakes with very different characteristic spatial scales along the flow and in 
transverse directions. In optics quite similar structures are called caustics [8]. The 
simplest way to represent such structures is to consider them as disks (singular 
manifolds) on which the density undergoes jumps vanishing at the disk boundary. 
It is worth noting that in classical hydrodynamics these regions are considered as 
"forbidden" where the solution can not be constructed anymore in a well defined 
way, that means as a single value solution. There are special procedures of con- 
struction of so called "shock type" solution, see the book of Whitham [9] for more 

^In the case of three- wave interacting waves of different branches the criterion of weak turbu- 
lence consists in smallness of the inverse nonhnear time defining by the kinetic equations in com- 
parison with the maximal growth rate of the corresponding decay instability for the monochro- 
matic wave. The latter in such system represents the inverse time of randomization (see, e.g. 
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detail. However, in collisionless plasma physics the detailed description is written 
in terms of velocity distribution function that satisfies to Vlasov equation in the 
phase space (pi, ri). For this function the regions where the hydrodynamic solution 
becomes multi-valued and poorly determined are the areas where the distribution 
function undergoes the transition from the one single peak distribution to the one 
having three peaks, however there is no any crucial break in the phase space. From 
the other hand the spatial derivatives of some characteristics have infinite gradients 
also, but this effect can be explicitly described in the frame of the same Vlasov 
equation for collisionless plasmas. Formation of simple waves making use of such 
a description was considered by Gurevich and Pitaevskii pTO]. 

In weak turbulence waves are assumed to be randomly distributed with a weak 
correlation between waves because of weak nonlinear interaction. The process 
of breaking is purely coherent and, respectively, the expanding caustics should 
be treated as coherent elementary entities of wave turbulence. We will call the 
turbulence where this effect becomes important a moderately strong turbulence if 
the density of caustics is still small enough to neglect their intersections that are 
supposed to be rare. The caustics can be described if their centers and orientations 
are determined. We will assume their distribution in space to be random. From 
the other hand, the jumps as the density singularities are known to result in a 
power-law tails in the short-wavelength part of the turbulence spectrum. This 
idea was first proposed by Phillips |TT] and allowed him to determine the water- 
wave turbulence spectrum in the presence of whitecaps, i.e., of the singularities 
on the fluid surface. Later the very same idea was developed by Saffman [12] 
to determine the isotropic spectrum of two-dimensional hydrodynamic turbulence 
at high Reynolds numbers (for details see [131 [II])- Kadomtsev and Petviashvili 
(KP) [15] suggested that the acoustic turbulence can be considered as a randomly 
distributed set of shocks. For the isotropic case they found the spectrum of the 
energy distribution, now known as the KP spectrum, 

E{k) ~ k-^. (3) 

The goal of our paper is to determine how the KP spectrum is modified in the 
presence of strong enough anisotropy in plasma. Such an anisotropy can be en- 
tirely inherent, either due to the mean magnetic field or because of the anisotropy 
of turbulence source/pumping. It is worth noting that in the weak-turbulence 
regime, when the wave dispersion of acoustic waves can be neglected, the angu- 
lar distribution of the spectrum repeats the anisotropy of the pumping because of 
three- wave resonant conditions, ujk = + ^"^^ k = ki + k2. In this case 
the additional anisotropy of the spectra can appear either due to weak dispersion 
[ISl E] or because of nonlinear renormalizations to the wave kinetic equation [17] . 
For weak MHD turbulence, in both low and high beta plasmas, the situation is 
familiar to the weak acoustic turbulence in isotropic media for the nonlinear inter- 
acting fast magnetoacoustic waves which dispersion is mainly defined by modulus 
of the wave vector. However, at /? -C 1 slow magnetoacoustic waves undergo strong 
anisotropy in both dispersion and nonlinear interaction due to the external mag- 
netic field. As it was shown in the paper [18] by one of the authors, it results in 
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appearance of anisotropic Kolmogorov-type spectra with power dependences rel- 
ative to both longitudinal and transverse projections of the wave vector. Such 
character of the spectra can be established also for the interaction of Alfvenic and 
slow magnetoacoustic waves at low beta plasma [1]. For (3^1 slow magnetoa- 
coustic waves have the same dispersion as the Alfvenic waves. This degeneracy 
also changes significantly the nonlinear interaction between waves and leads in the 
weak-turbulence regime to the spectrum of the Kraichnan-Iroshnikov type [19j (for 
details see pUlfT]). 

In the present paper, for moderately strong acoustic turbulence, we study how 
the angular ordering of shocks can change the angular structure of the spectrum. 
We show that for the strong anisotropy the spectrum has the jet-type behavior 
with power increasing along the jet with the same exponent as for the isotropic 
KP spectrum and the falling-off dependence in transverse direction ~ k]_^. The 
latter originates from the contribution from boundaries of caustics (disks). It is 
necessary to mention that for two-dimensional acoustic turbulence the phenomenon 
of the jet-type spectra generated by shocks was observed first time in the numerical 
experiments ^IJ. However, comparably small spatial resolution could not allow the 
authors [2T] to treat the fine structure of jets. 

The second objective of our work is to show how the presence of the KP spec- 
trum and its relative role can be evaluated making use of real experimental data 
( for instance, from spacecraft data for solar wind turbulence) in the presence of 
weak turbulence tails falling off more slowly than the KP spectrum. This question 
is very important because in the case when weak and strong turbulent components 
co-exist, the KP component of the spectrum can provide a quantitative measure 
of the relative role of coherent structures for moderate acoustic type turbulence. 
Hereafter we shall also show that the spectral index of the isotropic KP spectrum 
([3D as a function of the wave-vector does not depend on space dimension. It is 
worth noting that the Fourier transform of the density single point auto-correlation 
function dependence upon time shows the same power-law spectral index in time 
domain (in frequency) as in the spatial domain (spectrum in k) [31 

2 Isotropic spectra 

Let us consider the contribution of shocks to the frequency spectrum, i.e., the 
energy distribution dependence on frequency. To find the spectrum one should 
calculate the auto-correlation function for the density p(t) measured at some point 
Tq as a function of time: K{t) = {p{t + T)p{t)) where the angular brackets stand for 
time averaging of the (mass) density p{t) and then carry out the Fourier transform 



Here the density distribution is assumed to be homogeneous so that K^^ does not 
depend on fq. It is worth noting that defined in such a way, coincides, up to 



K{t), 



K 
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constant factor, with the energy spectrum in the frequency domain: 

Po 

In the weak-turbulence approximation, in an isotropic case E^^ can be expressed in 
terms of the spectrum e{k) by means of the relation 

E„ = e(k) (4) 

where k = u/cg. This formula is the result of integration of the energy spectral 
density in the k — u representation Eku = £{k) 5{ijJ — Uk) with respect to k. In the 
weak turbulence regime the spectral density E^^^ has the (5-function dependence 
upon frequency indicating that the wave ensemble is weakly nonlinear. In the 
presence of shocks, i.e. for the strong turbulence regime, such relations are not 
valid any more. 

Our aim now is to determine the contribution of shocks associated with the 
density jumps to the E,^ spectrum. To achieve this one should take into account 
that at the instant ti of jump passage through the measurement point fq the first 
derivative dpjdt is proportional to 5{t — ti), i.e., 

— = ^ Api6{t — ti) + regular terms. (5) 
ot ■ 

Assuming that density jumps Api and crossing times ti are random quantities one 
can calculate the contribution of these singularities in Eq. to the spectrum. The 
Fourier transform of the contribution associated with these terms can be written 
as follows: 

p» = ^i:Ap.e-'"'-- (6) 



Here, 



oo 



Pu. = I p(t)e*"*rft, p{t) = ^ I e-'-'p^dcu. 



To find the E^ spectrum, one should square the absolute value of and average 
the resulting expression. The averaging over crossing times ti yields 



E. = 7r^{\pJ') = 7r^,i^Py^ (7) 

where u = N/t is the jump crossings frequency, here N is the number of dis- 
continuities met during the averaging time r, and (Ap)^ is the average value of 
(Ap)2. 

The very same approach can be applied for finding the spatial spectrum of 
acoustic turbulence for the one dimensional case [D = 1) when instead of ([5]) we 
have 

dp ^ 

— = ^ Api6{x — Xi) + regular terms (8) 
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where the positions of jumps along Hence we get the following ID 

spectrum E{k): 



Here, ni is the number density of shocks per unit length, po is the mean density 
of medium (per unit length), and (Ap)^ is the mean-square of the density jump at 
the discontinuity. 

This calculation can be considered in isotropic three dimensional C3iSG clS the 
estimate of the correlation function along any chosen straight line. Now we shall 
make use of it to obtain the three-dimensional (3D) isotropic spectrum ([9]). To this 
end one should notice that the density correlation function for isotropic turbulence, 
4'i.yi) = {p{^i + Vij ^2, X3)p{xi, X2, X3)) , has the Fourier spectrum (with respect to 
only one single variable yi\) that coincides, to within a factor, with Eq. (Q: 

Ni 



27rA;2 



(Ap)2, (10) 



where A'^i in this case should be considered as the mean linear density of discon- 
tinuities. This correlation function (pk is related to the three-dimensional Fourier 
spectrum 

$(|k|) = J (j){r)e-'^-''dr 



by the following formula 



//•oo 
$(ik|)dk^ = TT ^s)d^' 



where k^ is the component of the wave vector k perpendicular to the a;-axis. This 
allows one to obtain differentiating this equality with respect to ki the following 
relation: 

^ ^ 27ik dk ■ 

Then, substituting Eq. (fTOl) . one can find the following expression for the spectrum 

The same approach is applicable in two-dimensional case. It is a little more difficult 
technically because it requires to solve the integral Abel equation (e.g. compare 
with [T2I [I3]). However, making calculations one can find that the spectrum de- 
pendence upon k has the same spectral index as in one- and three- dimensional 
cases: 

E2ik) oc k-"^. 

According to our knowledge spectrum IQ was first obtained by Burgers in [22] 
and, its generalization for mult i- dimensional situation was found (|TTi) by Kadomt- 
sev and Petviashvili [15]. In the next section we show how to generalize it taking 
into account the effect of anisotropy. 
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3 Anisotropic KP spectra 



The analysis of anisotropic situation consists in taking account of two geometric 
factors. One is connected with the anisotropy of the caustics orientations, another 
with the anisotropy due to the "emission" diagram of each caustic that can be 
considered as the source of the wave spectrum. For the sake of simphcity and 
without loss of generality we can consider that each single caustic has the form of 
a disk. This simplification is based on important property of these objects, namely, 
during their evolution most of time their relative size in the perpendicular direction 
is much larger than along the parallel axis. 

Let us consider a single caustic with radius Ri perpendicular to the x -axis and 
centered at the point fq = (xq, r^o)- To consider this effect equation should be 
replaced by the following one: 

dp 

— = Ap(|r_L — r_Lo|)(5(x — xq) + regular terms. (12) 
ox 

Here, Ap{rj_) is considered as a continuous cylindrically-symmetric function of r± 
that vanishes at the disk boundary r± = R {Ap{R) = 0) and remains zero outside 
the disk. 

Then, the Fourier transform of the terms corresponding to singular part of Eq. 
( fT2|) is given by the integral 

? r 27r? r^' 

Jr^KB. JO 

where k = {kx,ii±) and Jo{k±rj_) is the Bessel function. This is the contribution 
from one single singularity. The total contribution from all discontinuities can be 
found as the sum: 

Pk = -2vrz V- / Ap{r±) r±Jo{k^ar±)dr^. 

^ kilo "'0 

Here, iIq is the normal unit vector to the discontinuity a, are the disk center 
coordinates and k±a is the transverse projection of the wave vector k to the disc 
plane {k\^ = k"^ — (kn^)^). It is worth noting here that the anisotropic character- 
istics of the spectrum is related to the anisotropy of the distribution of unit vectors 

To find the spectrum of turbulence, one should average \pk\^ over all random 
variables. It is natural to assume that the coordinates r,^ of centers of caustics are 
distributed uniformly, the averaging over these variables results in 



\Pk 



2 = 47r2Ar 



/ Ap{r^)r^jQ{k^ri_)dr^ 
Jo 



1 rR 

ku Jo 



(13) 



Here, k^^ = kn, is the density of discontinuities and the angular brackets stand 
for the averaging over characteristic sizes R and angles. 
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Our major interest here is in evaluation of the short-wavelength asymptotic 
behavior of the spectrum found (IT^ . i.e., kR S> 1, where R is the characteristic 
disc radius. Thus, if k^^R ^ 1, then the expression inside the integral in Eq. f|T3|) 
represents rapidly oscillating function (due to the Bessel function Jo{k±r±)) and 
therefore the integral can be evaluated by means of the stationary phase method. 
Using the relation 

xJo{ax) = — — {xJi{ax)) , 
a ax 

and integrating by parts we have 

rR R d 

I 'Ap{r±)rj_Jo{k±rj_)drj_ = —- — / Ji{k±Rx) x—Ap{Rx)dx 
Jo k± Jo dx 

where the property Ap{R) = is used. Having in mind that the major input to 
the integral comes from the vicinity of the boundary at = i? we can use the 
asymptotic expression for the Bessel function at large k±R 




and obtain the following asymptotic estimate of the integral: 

d f2Rsm{k±R) d 



f /r Ti N 1/2 " A /T-. X r ZKsm[K±K) a . / X 

/ cos [kj_Rx) X ' —Ap[Rx)dx ~ —\ ^ ; — Ap[r±) 

Jo dx V TT k/ dr± 




r^=R 



where the main contribution to the integral comes explicitly from the boundary 
= R. 

Hence the resulting spectrum e(k) (still prior to angular averaging!) can be 
written as follows 



(14) 



It is worth noting that the stationary phase method used for the evaluation of the 
integral in Eq. f[T^ is applicable at almost all angles 6 {6 is the angle between 

the vectors k and n) but not in two narrow cones 6 < ^o = {kR^ and it — 6 < 
i^o- Inside these cones one can find that the integral can be considered to be 
independent on k (suggesting cos {k±Rx) ~ 1). The spectrum e(k) inside these 
cones is then given by 

Po kl 



where 



rR 

r = / Ap{r±)r±dr±. 
Jo 



From Eq. (HM one can see that the spectrum ei(k) contains three singularities, 
namely, at angles 9 close to 0, vr and 7r/2. For angles close to the cone 6 ~ (jcR^ 
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and 71 — 9 ^ (jtRj expression (|T^ matches Eq.( JT5l) . For angles close to 7r/2, in 
Eq. f|T^ an additional effect should be taken into account, namely, the bending 
of the caustics. If a is the characteristic value of bending, Eq. ( fT4l) is valid in the 
region \9 — 7r/2| > (ka)^^. 

Distributions (HM and ( fT5i) allow one to carry out the averaging over angular 
distribution of caustics normals, if their distribution is anisotropic, i.e., calculate 
the angular dependence of the energy spectrum -E'(k) = A;^e(k) in short enough 
scales (kR) > 1. 

First, let us check that in the isotropic case the KP spectrum (11 II) follows from 
the above calculations. In this case averaging over the angles corresponds to the 
integration of Eqs. f|T^ and f|T5l) with respect to 6. The integration of expression 
(fT^ near poles 9 = 0,71 gives 

E,{k)=47re i,(k)9d9 = ^^^^^4J-. (16) 
Jo pok^ R 

In the integration of Eq. f[T^ over angles, the main contribution to the spectrum 
comes from the angles close to 0, tt and 7r/2, where the spectrum f[T^ has singu- 
larities. For 9^0 (for 9 ^ tc), the integration is cut off at angles 9^ = { at 
9k = IT — "^o), and for — ^ 7r/2 it is cut off at angles \7t/2 ± 9k\ ~ {ka)~^. As 
the result of averaging of the expression (IT^ we get that the main contributions 
coming from the angles near the cone: 

The spectrum E{k) is given by the sum of (ITBi) and (ITTI) . 



Pok^ \ R^ 



that has the same dependence on k as the KP spectrum f fTT]) obtained in a little 
different way from conventional consideration. 

Our above performed analysis of spectral angular dependence upon the angle 
that is applicable first of all to one single caustic leads to the conclusion that the 
"emission diagram" of it represents narrow cones around 6' = 0, |, vr. The angular 

width of these cones is proportional to (JvRJ around 9 = 0,it and to (ka)^^ 
around 9 = 7i/2. In the case of anisotropic distribution of caustics normals the 
angular dependence of the spectrum will simply reproduce this angular dependence 
if the width A9 of this angular distribution is sufficiently larger than the width of 
these cones. This can be resumed as follows, if the angular distribution of caustics 
normals is wide enough so that > then after averaging of e(k) one should get 
the same dependence of the spectrum upon k at large enough k as for the isotropic 
KP spectrum, i.e. ~ k~^ with the degree of the anisotropy exactly the same as the 
distribution of caustics normals. From the other hand if the angular distribution of 
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caustics normals is sufficiently narrow , e.g., if all falling shock fronts are oriented 
almost unidirectionally (this can be caused, e.g., by pumping or initial/boundary 
conditions), then the spectrum will have a sharp peak in this direction. If the width 
A6' of angular distribution is narrower than -^o, i.e., if A6 < -^o, then the spectrum 
E{k, 6) up to the multiplier k"^ will repeat the distribution given by Eqs. f|T^ and 
(|T5|) . namely, the distribution in fc-space will have the form of the jet . In the cone 
9 < the spectrum has a maximum with fall off at large k in accordance with 
([T^ . i.e ~ k'^'^ like the KP spectrum flTTl) . At larger angles, 9 > the spectrum 

E{k, 9) will rapidly decrease in the transverse direction proportionally to fcj^. 

Note, however, that this asymptotic behavior is intermediate, because do = 
{kR)~^ decreases with increasing k. For this reason, when averaging over angles, 
singularities in f|T^ become essential for 6' — > 0, and, starting from certain k = k*, 
the spectrum will decrease as k~'^ with increasing k. The angular width of the 
spectrum will be of order A9 . 
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